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Abstract 

We study spontaneous symmetry breaking in a one-dimensional 
driven two-species stochastic cellular automaton with parallel sublat- 
tice update and open boundaries. The dynamics are symmetric with 
respect to interchange of particles. Starting from an empty initial 
lattice, the system enters a symmetry broken state after some time 
Ti through an amplification loop of initial fluctuations. It remains 
in the symmetry broken state for a time T2 through a traffic jam 
effect. Applying a simple martingale argument, we obtain rigorous 
asymptotic estimates for the expected times < T\ >oc LlnL and 
In < T2 > oc L, where L is the system size. The actual value of Ti 
depends strongly on the initial fluctuation in the amplification loop. 
Numerical simulations suggest that Tb is exponentially distributed 
with a mean that grows exponentially in system size. For the phase 
transition line we argue and confirm by simulations that the flipping 
time between sign changes of the difference of particle numbers ap- 
proaches an algebraic distribution as the system size tends to infinity. 

Keywords, spontaneous symmetry breaking, bridge model, cellular 
automaton, two-component exclusion process, martingale 



1 Introduction 

Spontaneous symmetry breaking (SSB) is associated with phase transitions 
and is therefore not expected at positive temperature in one-dimensional 
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equilibrium systems with short-range interactions. The underlying physical 
picture behind the absence of 1-d phase transitions, viz. the unsuppressed 
creation of islands of the minority phase inside a region of the majority 
phase (e.g. in an Ising system) due to thermal noise, is very robust. So 
it came as a bit of a surprise that in a quite simple stochastic (and hence 
noisy) lattice model of a driven diffusive system with short-range interac- 
tions SSB was observed Q3H]- In this so-called bridge model, two species 
of particles A, B move in opposite directions with rate 1 and are injected 
with rate a and ejected with rate (3 at the boundary sites. Although the 
dynamical rules are symmetric with respect to the two species, two phases 
with non-symmetrical steady states were found in a mean-field approxima- 
tion. In the symmetry broken states there is a macroscopic excess amount 
of one particle species, i.e., the order parameter A = pa~ Pb measuring the 
difference of the average particle densities of the two species A, B attains 
a non-zero value in the thermodynamic limit L — > oo. These analytical re- 
sults were confirmed by Monte Carlo simulations of finite systems of large 
size L. 

In the limit of vanishing ejection rate (/? — > 0) the existence of SSB in 
this model could be established rigorously [3J. It emerged that (at least in 
this limit) the phases of spontaneously broken symmetry are dynamically 
sustained by a traffic jam effect: The particles of one species pile up at one 
end of the chain (because of the small rate of ejection) and thus prevent 
the entrance of particles of the other species, until by an exponentially rare 
fluctuation (i.e. with a probability exponentially small in system size) no 
particles of that species enter for sufficiently long time. Then the traffic jam 
dissolves, allowing for particles of the other species to take over. Later some 
other stochastic 1-d lattice gas models exhibiting SSB were discovered [H[S], 
but the nature of the phase transition in the bridge model has remained 
obscure [H[7]. There is, in fact, recent numerical evidence suggesting that 
one of the two symmetry broken phases vanishes in the thermodynamic 
limit [E1[H]. 

It would seem natural to attack the problem of SSB from a macroscopic 
viewpoint by deriving a hydrodynamic description of the lattice gas model 
under Eulerian scaling. This approach indeed works for vanishing boundary 
rates [TO] , but fails for the general case due to the lack of a sufficiently 
general hydrodynamic theory for two-component systems in the presence 
of boundaries [llj . Only partial results for some specific infinite multi- 
component systems are known [T21 [J3J [JJl [J51 [TBI 113 ■ 

These and other puzzles make driven diffusive two-component systems 
a matter of considerable current interest, see [TS] for a review. In [Jj5] we 
studied a variation of the bridge model with parallel sublattice update. 
The deterministic bulk update scheme simplifies the treatment of parti- 
cle transport, while stochastic creation/annihilation events occur at the 
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boundaries in a similar fashion as in the original bridge model. Thus - 
while maintaining noisy dynamics - the complexity of the problem is re- 
duced. This allowed us to determine the exact phase diagram (for all values 
of the boundary rates) and to elucidate the mechanism that leads from a 
symmetric particle configuration into a state with broken symmetry. In this 
work we make rigorous the main results reported in [19j . which are com- 
plemented by new heuristic results on the transition line. We also point 
out that earlier work on a similar single-component system |20j yields a 
rigorous asymptotic estimate for the residence time in a symmetry broken 
quasi-stationary state of the two-component system and we present new 
results on the phase transition between the symmetric and the symmetry 
broken phase. 

In Sec. 2 we define the model and state our main results for the sym- 
metry broken phase. The proofs are given in Sec. 3. In Sec. 4 we present 
simulation data which provide further insight in the relevant time scales of 
the model. Analytical and simulation results for the phase transition line 
are given in Sec. 5. We conclude in Sec. 6 with some brief remarks. 

2 Model and results 

2.1 Bridge model with sublattice parallel update 

The model considered here is defined on a one-dimensional lattice of length 
L, where L is an even number. Sites are either empty or occupied by a 
single particle of either species A or B, i.e., the particles are subject to an 
exclusion interaction and the occupation numbers r]^(i) and ??b(«) of each 
site i obey Tjji(i) + < 1. The dynamics is defined as a parallel sub- 
lattice update scheme in two half steps. In the first half-step the following 
processes take place: At site 1 it is attempted to create a particle of species 
A with probability a S [0, 1] if the site is empty, or to annihilate a particle 
of species B with probability (3 £ [0, 1], provided the site is occupied by 
such a particle: 



Accordingly, at site L a particle of species B is created with probability a 
and a particle of species A is annihilated with probability (3: 



In the bulk, the following hopping processes occur deterministically be- 
tween sites 2i and 2i + 1 with < i < L/2: 




(1) 




(2) 



AO^OA , OB^BO , AB^BA . 



(3) 
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Figure 1: Stationary phase diagram of the sublattice bridge model as a 
function of creation and annihilation probabilities a and (3. 



In the second half-step, these deterministic hopping processes take place 
between sites 2i — 1 and 2i with < i < L/2. Note that the dynamics 
is symmetric with respect to interchange of the two particles species com- 
bined with space reflection. The original bridge model [TJ [5] arises as the 
continuous-time limit of this model with stochastic hopping and has the 
same symmetry. 

2.2 Results 

The stationary phase diagram of the model in terms of the parameters a 
and (3 is presented in Figure Q] The boundary lines are not difficult to 
analyze. Along the line a = 0, < /3 < 1 there is no injection and the 
stationary state is the empty lattice (which is trivially symmetric in A and 
B). For /9 = 1, 0<a<litis easy to verify by direct computation that the 
product measure with alternating densities Pa(}) — 0, Ps(i) = ol/{\ + a) if 
i is odd, and pa(i) = J r a )i PB(i) — if i is even is invariant. Also this 
stationary state is symmetric. For /3 = 0, < a < 1 there is no ejection and 
the system is highly non-ergodic. Any blocking measure with A-particles 
accumulating at the right boundary and -B-particles accumulating at the 
left boundary is invariant. We notice that most of these measures are not 
symmetric. Finally, for a = l,0</3<l there are two stationary product 
measures, one with = for all i and pa{i) = 1 for even i, pa(i) = 1—0 

for odd i and an analogous one with A and B particles interchanged. In 
each case one species is completely expelled from the system even on a 
finite lattice. This is a trivial, absorbing form of SSB with no transition 
between the two states. The dynamics reduce to that of the single-species 
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a = 0.6, /3 = 0.2 
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Figure 2: Average density profiles obtained from Monte Carlo simulations 
in the symmetric phase (top left), the broken phase (top right) and on the 
transition line (bottom). A densities are shown by o, B densities by ■. 



lattice gas studied by one of us earlier [2D]. In the special deterministic 
point a — (3 = I there are three invariant measures arising as limits of the 
measures described above. One is the symmetric alternating A/B measure 
{13 = 1, a — > 1), the other two are the symmetry-broken measures with only 
A or only B particles (a = — > 1). These considerations hold for all 
system sizes L, and in what follows the boundary lines are excluded from 
the discussion. 

The interior of the phase diagram can be explored by Monte Carlo 
simulations [TDJ. Two phases are found: 

• If a < (3, the system exhibits a symmetric steady state. Here, the 
bulk densities in the limit L — > oo are pA(i) = 0, ps(i) = a{3/(a + j3) 
if i is odd, and pa(}) = a/3 /(a + f3), ps(«) = if i is even. 

• If a > [3, the system resides in the symmetry broken phase. Assume 
the A particles to be in the majority Then, the bulk densities in the 
limit L — > oo are ps(i) = for all i, pA{i) — 1 f° r i even and pa{^) 

1 i for i odd. This means that the symmetry is maximally broken 
and the minority species is completely expelled from the system. 

• On the transition line for a = (3 the system switches between sym- 
metric and broken states, which is studied numerically in Section 5. 
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The intention of this paper is to further elucidate the dynamics leading to 
symmetry breaking which were identified in [19] and to determine rigorously 
the time scales associated with SSB for a > (3. 

If one of the species is expelled from the system, the dynamics of the 
majority species for a > [3 is identical to the single species ASEP with 
parallel sublattice update. The exact steady state distribution fi for this 
system has been characterized in [20] for all a, (3 and all system sizes L. 
Let Sq be the distribution that concentrates on the empty lattice for the 
single species system. Then we write 

v A = M ® an d v B = 8q® \x (4) 

for the distributions of our two species system, where one of the species 
is expelled and the other distributed according to fi. Note that we do 
not explicitly write the dependence of va and vb on a, f3 and L. We 
say that the system reaches a configuration r\ — (77^(1), rjB(i)) i _ 1 L with 
broken symmetry if rjA(i) = 1 for all even i and r)B(i) = for all i, or 
vice versa t)b{i) = 1 for all odd i and rjA{i) = for all i. Let Ti(tj) 
be the (random) time when the system first reaches a symmetry broken 
configuration, starting from configuration r). T2(rj) is defined to be the 
time until the first particle of species B enters the system. With these 
definitions we can state the main results of this paper. 

Theorem 1 Time to reach a symmetry broken configuration 

For a > (3 we have lim sup — ^— — < 00 . 

l^co LmL 

Here (..)„ denotes the expected value over the time evolution of the 
process with initial distribution v. So if we start the process with an empty 
lattice (So = Sq ® Sq) then the expected time to reach a symmetry broken 
state is not growing faster than LlnL with the system size L. 

Theorem 2 Stability of states with broken symmetry 

For a > (3 we have lim inf ^ n ^ 2 ^ VA > q _ 

L — too L 

So starting with a symmetry broken state va where the B particles 
are expelled from the system, the expected time until the next B particle 
enters grows exponentially with L. By symmetry an analogous statement 
holds if A and B particles are interchanged. 

Since va = A ® ^0 ana - I 1 1Sl stationary for the single species system, 
va is invariant for the two species system for times t < T2. So although 
va and vb are not stationary for finite L, they are exponentially stable 
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by Theorem 2 for a > (3. Together with Theorem 1 this provides a proof 
for spontaneous symmetry breaking in this model in the interior region 
of the phase diagram. Note that this argument is done without knowing 
the stationary distribution exactly, which should of course be close to the 
mixture ^A+h^B, corresponding to the usual concept of a phase transition 
in the context of Gibbs measures (e.g. in an Ising system). 

The proofs are given in Section 3, where the proof of Theorem 2 is a 
straightforward application of the results in [20]. The proof of Theorem 
1 relies on a simple martingale argument for an interesting amplification 
mechanism which has been published in [19] . and will be explained again 
in the next subsection for self-containedness. 

2.3 Dynamics of symmetry breaking 

It is assumed that at t = Q there are no particles in the system and that 
a > f3 > 0. Other symmetric initial conditions with a non-empty lattice 
can be treated in a similar fashion. Starting from the empty lattice, A 
(B) particles are created at every time step with probability a at site 1 
(L). Once injected, particles move deterministically with velocity 2 (—2). 
Therefore, at time t = L/2 the system is in a state where the density of A 
(B) particles is a (0) at all even sites and (a) at all odd sites. In this 
situation both creation and annihilation of particles are possible. 

However, it turns out that the effect of creation of particles is negligible: 
Since a > (3 the deterministic hopping with velocity 2 transports on average 
more ^4-particles towards site L than can be annihilated there. This leads 
to the formation of an A-particle jam at the right boundary, blocking 
the injection of i?-particles. An analogous argument holds for the left 
boundary, which is blocked by a B-particle jam. In these jams, the only 
source of vacancies is annihilation at the boundaries with probability (3 
in the first half-step. In the second half-step the vacancy moves one site 
towards the bulk with probability 1. Therefore, in a jam, the density of A 
(B) particles at even (odd) sites is 1, while that at odd (even) sites is 1 — (3. 
So the only way to create particles in this situation is a complete dissolution 
of a jam. But as long as it gains particles from transport through the bulk 
this is a very rare event since a > f3. We show in Lemma [3. 21 (Section 3) that 
the average number of created particles is small and bounded independent 
of L. So creation of particles in this jammed situation becomes negligible 
in the limit L — > oo and will be neglected in the following explanation. 

The number of particles in each of the two jams reduces by one in 
every time step with probability (3. Since creation of particles is negligible, 
the influx into the jam ceases after some time and the jam eventually 
dissolves. By fluctuations, one of the jams, say the B-jam at the left 
boundary, dissolves first. A particles can enter the system while B particles 
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are still blocked until the A-jam at the right boundary is also dissolved. 
The configuration of the system at this point is illustrated in Figure [3] at 
time t\ (setting k = 1). The light grey region denotes a region of low 
density of A particles where the density is a (0) on even (odd) sites. The 
(random) length of this region, Ai±, describes the majority of one of the 
species. The description is symmetric, so if A£± < this corresponds to a 
majority of B particles. Thus the average value (A£i)s = 0, but typically 
A£\ = 0(yL) due to fluctuations for large L and one of the species has 
the majority (see Lemma 13.31 in Section 3). 

The time evolution just described constitutes the first cycle of a peri- 
odic behavior which can be effectively described by the dynamics of low 
density regions and jams at the boundaries. The key ingredient for this 
simplification is the jamming mechanism described above. The cyclic be- 
havior consists of 4 stages, which we summarize in the following and which 
are illustrated in Figure [3l For simplicity of presentation we assume that 
A£fe > 0, i.e. the A particles have the majority. 

1. At the beginning of a cycle (t\ := 0) there is a low density region of 
A particles at the left boundary of length \A£k\ > 0. 

Both species enter the system with probability a and penetrate the 
bulk dcterministically with speed 2 (—2). 

2. The low density region of A particles reaches the right boundary at 
time t2 := (L — |A£fc|) /2, blocking the creation of B particles. 

A particles still enter with probability a and exit with probability 
/3 < a, further increasing their majority. 

3. At time ts :— L/2 the B particles reach the left boundary, blocking 
also the creation of A particles. 

Both species form jams at the boundaries, which gain particles from 
the low density regions. Since creation of particles at the boundaries 
is negligible, both jams eventually dissolve. 

4. Let £4 be the time when the jam of B particles is dissolved. 

A particles start to enter the system. Again, since a > (5 the majority 
of A particles increases on average. 

5. At time t§ the A-jam at the right boundary is dissolved and also B 
particles can enter the system. 

The cycle is finished when both jams are dissolved, i.e. at time max{£4, t$}. 
Note that for given Ai^, t\, i2 and £3 are deterministic, whereas £4 and t$ 
are random times, which will be defined more precisely in Section 3. In 
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Figure 3: Cyclic behavior of the dynamics of symmetry breaking. 
Top: Illustration of the stages involved in the fc-th cycle of the amplification 
loop as explained in the text. Low density regions are drawn light grey, 
jams are dark grey and white regions of the system are empty. On the 
right the stages are identified in a blow-up of simulation data shown in the 
bottom. 

Bottom: Two realizations of MC simulation of symmetry breaking, starting 
from the empty lattice with a = 0.9, (i = 0.8 and L = 10000. The density 

p^is drawn in a full line ( — ), ps in dashed line ( ) and the dotted line 

indicates the density in the symmetry broken state. 
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Figure [3] it is assumed that the B-jam dissolves first, i.e. t± < t 5 , which is 
most likely if A£ k > 0. But in general £4 > t$ is also possible and included 
in the above description. The result of the cycle is 

A4+1 := 2 (i 5 - U) , (5) 

which is the initial condition for the next cycle. As long as — L < A£ k +i < 
L the cycle can restart with < <2 < £3, and all the stages are well defined. 
Note that within this framework the process starting from the empty lattice 
is a cycle with initial condition A£ = and t% — t%. The cyclic behavior 
leads to an amplification of the initial fluctuation A£i , namely 

(A£ k \A£ 1 ) So ~ (2a/j3 - l) k Ah for large L , (6) 

as given in Lemma \'S. 2 1 Together with the fact that the average duration 
of a cycle is of order L (see Lemma [3.3[) . this is the core of the proof of 
Theorem [T] provided in the next subsection. 

3 Proofs 

3.1 Proof of Theorem 1 

In order to define the random times £4 and t§ in the cycle described above 
more precisely, we use the following procedure: At time £2 we mark the last 
particle of the minority species, corresponding to the rightmost B particle 
in the above description. If due to a fluctuation there is no particle of 
that species in the system, the cycle is finished at £2 with A£ k +\ = L and 
also the amplification loop stops. Otherwise, we analogously mark the last 
particle of the majority species at time t% and define t± and t§ to be the 
time when the marked B and A particle, respectively, just left the system. 
If there were no particles of the majority species in the system at time £3, 
we set £5 := £3 (or t± := £3 if the B particles have the majority). At time 
max{i4,is} we restart the cycle with initial condition A£ k +i given in ([S]). 
This determines a process (A£k)k=o,i,.. 011 Z with A£ = 0. If |A^ + i| > L 
the process stops and we reach a symmetry broken state. 

We first note some Lemmas used in the proof. The proofs of the Lemmas 
are given in the next subsection. 

Lemma 3.1 If \A£k\ > L the system reaches a symmetry broken configu- 
ration as defined in Section 2.2 within a time of order L. 

We define k to be the (random) number of cycles performed when the 
loop stops, i.e. 

k := inf {fc : \A£ k \ > L) . (7) 
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This is a stopping time for the process (A4)fe=o,i,.. and in the following 
we aim to find an estimate on the expected value (k). We will use the 
following recursion relation: 

Lemma 3.2 For a single cycle 

(A4+i | A4> = qM k + C k , (8) 

where q = (2 ^ — 1J > 1 and the constants are bounded independent of 
L and k. This constitutes an amplification of the initial fluctuations A4 
and we get inductively 

k 

(A4+i | Mi) = M iq k + c *V k ~ % • ( 9 ) 

Here and in the following we omit the subscript do since all expectations 
are understood with respect to the empty initial condition. 

Lemma 3.3 The average length of a cycle is bounded above by 

(i 5 |A4) = Q + l) \ + Q - l) <A4 I A4>/2 = 0(L) . (10) 

and the initial fluctuation is of order y/L since 

VariMr) = 4L^p® . (11) 

To get an estimate for (fc) we use the optional stopping theorem for mar- 
tingales. Since q > 1 and A4 = 0(y/~L) as shown in (jlip . it is clear from 
([U that (A4)fc is a sub-martingale for sufficiently large L. Conditioned on 
the initial value A4 we define the process (Yfc)fc=i,2,.. by Y\ := A£i and 

Y k - F fc _i := A4 - (A4 | A4-i) for k > 2 , (12) 

following Doob's decomposition of sub-martingales. By construction, (Yfe)fc 
is a martingale and k is a stopping time for (Yfc)fc. Thus by optional 
stopping and using Lemma 13.21 we have for k > 2 

= (y s |A4) - <n_ 1 |A4) = 

= (A4|A4)-9<A4_ 1 |A4)-<C S _ 1 |A4) • (13) 

Also by Lemma l3~2l we know that | (A4I A4) | _• qL + 0(l), since otherwise, 
the process would have stopped before k. Using © and (Cj,_ 1 \Ali) = 0(1) 
we get 

gi + 0(l)>g|(A4_ 1 |A4)h(|A4|+0(l))(<z S - 1 |A4) . (14) 
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Since q > 1 Jensen's inequality yields 

^'• >< - ^mrsk ■ (15) 

This leads to 

lnfi/|A4| +o(l)) 
fcA4 <2+ . 16 

In q 

Again with Jensen's inequality we have 

(\nLf\Ah\) <ln(i(l/|A4|)) (17) 
Since with LemmaGEl (l/|A4|) = 0(£~ 1/2 ) we get 

(fc)<0(l) + ^(l + o(l)) • (18) 

Therefore, taking the expected value w.r.t. A^, the expected total time 
spent in the amplification loop is of order LhiL. Together with Lemma 
13. II this finishes the proof. □ 



3.2 Proof of Lemmas 

Proof of Lemma 13.21 In the following we analyze the distributions of 
the random variables and t§ to get the time evolution of At^. Let t„ be 
the (random) time it takes for a jam of n particles to dissolve. With this 

U = L/2 + t Nb + E% , t 5 = (L- A4) /2 + t Na + E\ . (19) 

Here Na (Nb) denotes the number of A (B) particles that entered the 
system up to time £3 (£2) before blocking, including Alk- Ea (Eb) denotes 
the number of time steps where the A (B) jam is dissolved, i.e. site L (1) is 
empty before £4 (is), when the respective marked particle exits. These are 
fluctuations and may lead to single B (A) particles that enter the system 
due to lack of blockage. We call such particles discrepancies and below we 
show that their expected number is bounded independent of the system 
size L. Apart from that, the boundary site in a jam is always occupied and 
particles are annihilated with probability j3. So the time r € {1, 2, . . .} for 
one particle to leave the jam is a Geo{f3) geometric random variable with 

P(t = k) = 13(1 - (3) k - 1 , <r) = 1/0 . (20) 

Let t 1 for i = 1, . . . , n be n independent copies of r and n £ N an indepen- 
dent random variable. Then for 

n 

t„ = J2 t * we have < r "> = <">( r > = ^if 3 ■ ( 21 ) 

i=l 
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In p9[) n = Na or Nb is a i?i(i, a) binomial random variable with 



P(n = k) = Qa'fl-a)" , fe € {0, ...,£} , (n) = at . (22) 

For n = iV^ it is t = A4/2 + t 3 - B A and for n = iV B , t = t 2 - B B , 
where Ba (Bb) is the number of time steps where the entrance of A (B) 
is blocked by singular B (A) discrepancies. With (|2"2"|) we have 

(N A ) = (A4/2 + 1 3 - (B k A )) a = (L + Al k ) a/2 - a(B k A ) , 
(N B ) = (h - (B%)) a=(L- A£ k ) a/2 - a(B k B ) , (23) 



where all expected values are conditioned on A£ k - According to (I20|) . di- 
viding ([2^)1 by (3 yields (tn a ) and (tn b ). Using this and (f2"3")) the average 
value of A£ k +i = 2(t 5 — i 4 ) conditioned on A£ k is given by 

(A4+i | A4) = (2^ - l) A4 + 2(^ - £*) - 2^(B k A - B|) . (24) 

This is equal to §8§ with 

C k = 2(E A -E B )-2~(B<A-B k ) = 2(E k -E k )-2^(E A - l -E B - 1 ) (25) 

since the blocking of A particles is caused by discrepancies of the previous 
cycle, i.e. (B%) = ^(E^ 1 ) and analo gous for B particles. ([9]) follows 
directly by induction. 

In order to finish the proof it suffices to show that (E\) and (E B ) are 
bounded independent of L for all k. Recall that an A-jam is defined as a 
region of A particles at the right boundary with densities 1 on even and 
1 — /3 on odd sides. Denote by M t , t > t%, the number of particles in 
the j4-jam at time t. If Mt > it decreases by 1 with probability in 
each time step. Until the marked A particle reaches the jam, say at time 
t*(L), Mt increases at least by 1 with probability a > f3 in each time step. 
Due to the sublattice parallel update, also an increase by two particles 
is possible. Both statements are true only modulo finite corrections due 
to discrepancies. Nevertheless, M t performs a biased random walk and is 
increasing on average. Thus it visits only finitely often for t — > 00, and 
thus also for t < t*(L). Now 

E\ = \{t€ [t 2 , **(£)] : M t = 0}| , (26) 

and thus (E A ) is bounded independent of L. The same argument holds for 
(E B ), finishing the proof. □ 
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Proof of Lemma EOl With (19]) and ([23} we have 

(t 5 ) = (L - (A4»/2 + (L + (A4)) a/(2/3) - a<B£>//3 + , (27) 
where all the expected values are conditioned on A£\. This leads directly 

to ms]). 

In the initial cycle starting with the empty lattice, t\ and £5 are i.i.d.r.v's 
and thus with A£\ = 2(t& — £4) we get 

yar(A4) =4V r ar(i 5 -i 4 ) =8yar(t 4 ) • (28) 

Analogous to (f2"3"|) this is easily computed as 

Var{U) = (N A ) Var{T l ) + (t 1 ) 2 Var{N A ) , (29) 

where Na and Tj are defined as in the proof of Lemma 13.11 Thus with 
(N A ) = aL/2, Var(N A ) = a(l - a)L/2 and (t 4 ) = 1/0, Var^) = (1 - 
/?)//? we get 

7ar(A£ 1 )=4L ^ 2 -^- /3) (30) 

□ 

Proof of Lemma 13.11 Let A£k > L. Then the low density region of 
A particles extends over the whole lattice, and there is only a finite num- 
ber of B particles (discrepancies) in the system. Therefore an A-jam will 
form at the right boundary blocking the entrance of B particles, whereas 
A particles will not be blocked. The number of particles in the ^4-jam 
performs a biased random walk increasing on average as explained in the 
proof of Lemma 13. 2\ and thus the jam will reach the left boundary in an 
expected time of order L. At this point the density of A particles on all 
even sides is 1 and if there are any singular B particles left, they will leave 
the system in a time of order L and the system reaches a symmetry broken 
configuration as defined in Section 2.2. □ 

3.3 Proof of Theorem 2 

Assume the system to have symmetry broken distribution va = M ® So as 
defined in Section 2.2 with particle species B expelled from the system. 
The density profile of A particles is given by the stationary measure [i for 
the single species system which is known exactly [20 . For L — > 00 this 
means pA{i) = 1 + o(l) for even sites and Pa{i) — 1 — f3 + o(l) for odd 
sites, up to boundary effects at the left boundary with i = 0(1). Species 
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B is expelled from the system and injection of B particles is only possible 
if site L is empty. Exact expressions given in [20] (equation (18)) yield 



Thus the probability that site L is empty is exponentially small in the 
system size, and the expected time {T2) VA until the minority species can 
penetrate the system started in the broken state is exponentially large in 
L. This is not surprising even without knowledge of the exact expressions, 
since for injection of the first B particle the complete jam of A particles 
has to be dissolved against the drive a > (3 and this jam consists of the 
order of L particles. 



4 Simulation results 

4.1 Results for T x 

In the proof of Theorem Q] we identified a bound on the expected number 
of cycles (k) until the amplification loops stops, which is called k* in the 
following. This bound grows like k*(L) ~ C + hxL/{2\sxq) with increas- 
ing L (p~8|) . To estimate the constant C we replace (l/|A^i|) in (fTBj) by 
l/y/Var(Ah) = 4L a(2 ~" -/3) as given in JTIJ and get 

U*(T\ lnL i o i M/? 2 /(2a(2-q-/?)) (v> , 
2 In q zmq 

With this choice of C, k*(L) is no longer a strict upper bound but in very 
good agreement with simulation results for (fc), as can be seen in Figure 
[4] In particular the simulation data show the same logarithmic growth 
with prefactor (21ng) _1 as k*(L), so in this sense the rigorous bound of 
Theorem [T] is sharp. 

As can be seen in Figure [3] (bottom) the number of cycles needed in 
each realization depends heavily on the initial fluctuation A£±, which is 
a random quantity even in the limit L — * oo. Thus we do not expect a 
law of large numbers for k, which is also supported by Figure [4j The bars 
indicate the standard deviation of the distribution of k which is more or 
less independent of the system size. 

4.2 Stationary results 

The stationary dynamics consists mainly of flipping between the two sym- 
metry broken states which are close to va and vb- Wc define the flip times 
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Figure 4: Simulation results for (k) for different system sizes L. The data 
are obtained as averages over 100 realizations, the average values are shown 
as for a = 0.5, ft = 0.4 and ♦ for a = 0.9, j3 = 0.5. Errors are of the size of 
the symbols and the bars denote the standard deviations of the distribution 
of k. The full lines ( — ) give a linear fit to the data points which agree very 
well with the estimate k*(L) given in ([52]) . shown as dashed lines ( ). 




r to be the times between two consecutive sign changes of pa ^ Pb, where 
PA = \ Yli=i VA(i) is the (time dependent) number of A particles in the 
system normalized by L, and ps is defined analogously. A typical trajec- 
tory of (pa — ps)(t) is shown in Figure [5] (upper left), showing that there is 
a clear timescale of flipping between the two states. The lower left plot of 
Figure [5] shows simulation data for the cumulative tail of the distribution of 
the random variable r. This distribution clearly consists of two parts, one 
of which are small flip times r = 0(1) which result from fluctuations during 
a single transition between the two symmetry broken states. The relevant 
flip times are the ones that increase with the system size marking the life 
time of the symmetry broken states. In the upper right plot of Figure [5] 
the average value of these relevant flip times, denoted by (r)' is shown to 
increase exponentially in the system size as (r)' ~ z L . For a = 0.5, = 0.4 
we measure z = 1.31 ±0.02 and for a = 0.9, /3 = 0.5 z = 2.55 ±0.08. Both 
values are larger than a//3, respectively, consistent with the result (I31| for 
T2, which is a lower bound for the relevant flip times. Normalizing the 
data by (t)' in the lower two plots of Figure [5] results in a data collapse 
for the extensive part of the distribution. The non-extensive part collapses 
to a jump at x = in the cumulative tail (lower left plot), showing that 
a substantial fraction of sign changes have flip times of 0(1). The lower 
right shows a logarithmic plot of the tail of the renormalized extensive part 
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Figure 5: Study of the flip time distribution. The upper left plot shows 
PA — Pb as a function of time for a = 0.5 and (3 — 0.4 and L = 50. The 
other plots show data for a — 0.5, (3 — 0.4 (unfilled symbols) and a = 0.9, 
= 0.5 (black symbols). The average values of the relevant flip times (r)' 
(see text) are shown to be exponentially increasing in L in the upper right 
plot. The full lines give a linear fit to the data points. The lower left shows 
the cumulative tail of the distribution of normalized flip times t/(t)' for 
L = 30 (A), L = 40 (0), L = 50 (□) and L = 14 (A), L — 16 (♦), L = 18 
(■). The same symbols apply in the lower right, a logarithmic plot of the 
cumulative tail of the renormalized distribution of relevant flip times. The 
full line denotes an exponential distribution with parameter 1. 



of the distribution, i.e. the distribution of relevant flip times (denoted by 
P') showing good agreement with an exponential distribution of parameter 
1. Therefore we conclude that the relevant flip times have an exponential 
distribution, with average value increasing exponentially in L. 

5 Dynamics on the transition line 

For the borderline case a = f3 the dynamics of the system can still be effec- 
tively described in terms of boundary jams and low density regions. The 
cyclic behavior can be observed (see Figure El upper left), but fluctuations 
are larger since the end of a jam is diffusing and the cycle lengths, though 
still of order L, are strongly fluctuating. But according to (|24|) there is no 
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amplification of fluctuations during a cycle. Instead, A£k is not driven to- 
wards ±L but is diffusing, so a symmetry broken state can still be reached 
within 0(L 2 ) cycles, and thus T\ = 0(L 3 ). On the other hand, when the 
system is in one of the symmetry broken states, the length of the jam of 
the majority species is only diffusing. So it dissolves in a time of only 
T2 = 0(L 2 ), which is the lifetime of a symmetry broken state for a = f3. 
Thus, no symmetry breaking takes place in this case. Instead, for large L a 
typical configuration is taken from a cycle, consisting of jams with diffusing 
length and of low density regions for both species. An average over many 
realizations leads to an approximately linear stationary density profile as 
shown in Figure [1] (lower left). Further, for a — (3 site 2 (L) is occupied by 
A particles for approximately half of a cycle length with probability a (1), 
leading to the stationary densities pa{2) = a/2 and pa{L) = 1/2. For odd 
sides an analogous argument yields pa(X) = and pa(L— 1) = (1— a)/2, 
which agrees well with Figure [T] 

Moreover, the formation and dissolution of boundary jams for the two 
species shows interesting temporal correlations |llj . In the following we 
study the distribution of flip times r between sign changes of {pa — 
analogous to the symmetry broken case. In Figure [HI we show three plots 
of pa — pb against time on different time scales for a = (3 = 0.5. For these 
parameters the maximal difference, given by the average density in the 
symmetry broken phase, is given by 1 — (3/2 = 0.75. In a time of order L 2 
the path explores the whole interval [—0.75, 0.75] (see lower right plot), but 
for smaller time scales the paths show a self similar structure (see lower left 
and upper right plot). Therefore one observes flip times on all length scales 
and in the limit L — > 00 we expect a scale free distribution of t. This is 
confirmed in Figure [71 showing a double logarithmic plot of the cumulative 
tail P(t > x). Data for different system sizes collapse without scaling and 
show agreement with a power law tail with exponent —0.5. For large x the 
data deviate from this behavior, due to boundary effects for smaller values 
of L as discussed above, and due to numerical inacurracies, using quantilcs 
to determine the cumulative tail. 

The exponent of the power law can be predicted by comparison with a 
random walk. L(pa — Pb), the difference in the number of particles as a 
function of time performs a symmetric random walk for a = (3 on the inter- 
val [— L,L]. Thus we use the scaling ansatz P(t > x) ~ L 1 g{x/L 2 ). The 
argument x/L 2 follows from the scaling of first passage times of symmetric 
random walks in one dimension and the fact that r is the return time to 
the origin. The power 7 can be fixed by the requirement that (r) = O(L), 
since the average return time to is inversely proportional to the stationary 
distribution at 0. This scales as 1/L because the stationary distribution of 
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Figure 6: Density profiles for a = (3 = 0.5 on different time scales. The 

upper left plot shows pa (full line — ) and pb (dashed line ). The other 

plots show the difference pa — Pb- The plots are stationary samples and 
do not start at t = 0, the axes only indicate the time scale. 
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Figure 7: Double logarithmic plot of the cumulative tail of the distribution 
of flip times r. Data for a = [3 = 0.9 are shown as ▲ (L — 3200), ♦ (L = 
6400), ■ (L = 12800), corresponding unfilled symbols for a — (3 = 0.5. 
The full line indicates the inclination of a power law exponent —0.5. 



19 



the random walker is a linear function on [-L, L}. We get 



(r)~|°°P^/i 2 )^Or 2 ), 



(33) 



and thus 7 = — 1 and a consistent ansatz is 



P(t >x)~ L^gix/L 2 ) . 



(34) 



To cancel the L-dependence for L — > 00 we have g(y) ~ as y — » for 

the scaling function and thus for x large enough 



giving the observed exponent —0.5. 

6 Conclusion 

In the present article we investigated spontaneous symmetry breaking (SSB) 
in a two-species driven cellular automaton model with deterministic bulk 
behavior and stochastic open boundary conditions. We analyzed in detail 
the dynamical mechanisms leading to SSB. Its main feature is a cyclic am- 
plification of initial fluctuations taking a time of order L In L, while a traffic 
jam phenomenon keeps the system in a SSB state for an exponentially large 
time. This lead to a proof of SSB in the thermodynamic limit using a sim- 
ple martingale argument without further assumptions on the rates, and to 
rigorous asymptotic estimates for the relevant time scales in the broken 
phase. The above mechanism is very different from the freezing-by-cooling 
scenario for broken ergodicity in one-component systems [4] that results 
from a localization of shocks [21j [22l [23] . 

Some dynamical and stationary properties at the phase transition line 
have been predicted analytically (but not rigorously) in terms of boundary 
jams and low density regions using the picture developed for the discussion 
of the broken phase. In particular, we found an asymptotically scale-free 
distribution of flip times between sign changes in the difference of particle 
numbers. The decay exponent of the distribution has been predicted us- 
ing random walk arguments and confirmed by numerical simulations. The 
exact phase transition line can be predicted correctly by a mean-field ap- 
proximation. The density profiles predicted in this way, however, differ 
from the numerically computed density profiles [19] . Therefore mean field 
theory is unreliable at the phase transition line, and was not presented 
in this paper, whereas the correct density profiles could be predicted in 
Section 5 by a rather simple argument. 




for x « L 2 , 



(35) 



20 



The amplification mechanism outlined above does not apply in the sym- 
metric phase a < (3 since the formation of boundary jams, a key ingredient 
for the amplification, does not work. The length of a boundary jam is 
driven towards small values so the boundary sites are not blocked and par- 
ticles are injected all the time. It is an intriguing open question whether 
similar mechanisms are at work in the original bridge model of [TJ [5] . 
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